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Abstract 

Methods in Riemann-Finsler geometry are applied to investigate 
bi-Hamiltonian structures and related mKdV hierarchies of soliton 
equations derived geometrically from regular Lagrangians and flows of 
non-stretching curves in tangent bundles. The total space geometry 
and nonholonomic flows of curves are defined by Lagrangian semis- 
prays inducing canonical nonlinear connections (N-connections), Sa- 
saki type metrics and linear connections. The simplest examples of 
such geometries are given by tangent bundles on Riemannian sym- 
metric spaces G/SO{n) provided with an N-connection structure and 
an adapted metric, for which we elaborate a complete classification, 
and by generalized Lagrange spaces with constant Hessian. In this 
approach, bi-Hamiltonian structures are derived for geometric me- 
chanical models and (pseudo) Riemannian metrics in gravity. The re- 
sults yield horizontal/ vertical pairs of vector sine-Gordon equations 
and vector mKdV equations, with the corresponding geometric curve 
flows in the hierarchies described in an explicit form by nonholonomic 
wave maps and mKdV analogs of nonholonomic Schrodinger maps on 
a tangent bundle. 
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1 Introduction 

Some interesting studies in geometric mechanics and field theory and the- 
ory of partial differential equations are related to nonholonomic structures 
characterizing integrability of certain nonlinear physical systems, their global 
and local symmetries [D [21 El IH E] ■ In parallel, the differential geometry of 
plane and space curves has received considerable attention in the theory of 
integrable nonlinear partial differential equations and applications to modern 
physics [6l [71 [HI m [TOl [Til [12] • More particularly, it is well known that both 
the modified Korteweg-de Vries (mKdV) equation and the sine-Gordon (SG) 
equation can be encoded as flows of the curvature invariant of plane curves in 
Euclidean plane geometry. Similarly, curve flows in Riemannian manifolds of 
constant curvature give rise to a vector generalization of the mKdV equation 
and encode its bi-Hamiltonian structure in a natural geometric way [131 IE] • 
This approach provides an elegant geometric origin for previous results on 
the Hamiltonian structure of multi-component mKdV equations (see, e.g. 
Ref. 

In recent work [17] , geometric flows of curves were studied in Riemannian 
symmetric spaces M = G/SO{n) which provide the simplest generalization 
of n-dimensional constant-curvature Riemannian geometries. The isome- 
try groups G of these spaces are exhausted by the Lie groups SO{n + 1) 
and SU{n), as known from Cartan's classification fT9]. The main results 
of [IT] were to show that, firstly, the Cartan structure equations for torsion 
and curvature of a moving parallel frame and its associated frame connec- 
tion 1-form encode 0{n — l)-invariant bi-Hamiltonian operators. Secondly, 
this bi-Hamiltonian structure generates a hierarchy of integrable flows of 
curves in which the frame components of the principal normal along the 
curve (analogous to curvature invariants) satisfy 0{n — l)-invariant multi- 
component soliton equations that include vector mKdV equations and vector 
sine-Gordon equations. The two groups G = SO{n + 1), SU{n) give differ- 
ent soliton hierarchies and account precisely for the two known integrable 
versions [IHl [201 [21] of vector mKdV equations and vector sine-Gordon equa- 
tions. Thirdly, the curve flows corresponding to such vector soliton equations 
were shown to be described geometrically by wave maps and mKdV analogs 
of Schrodinger maps into the curved manifolds M = G/SO{n). 
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A crucial condition [T7] behind such constructions is the fact that the 
frame curvature matrix is constant on these spaces M = G/SO{n). This ap- 
proach can be developed into a geometric formalism for mapping arbitrary 
(semi) Riemannian metrics [22] and regular Lagrange mechanical systems 
into bi-Hamiltonian structures and related solitonic equations following cer- 
tain methods elaborated in the geometry of generalized Finsler and Lagrange 
spaces [21 El [23] and nonholonomic manifolds with applications in modern 
gravity [H EHl EE] . 

The first aim of this paper is to prove that solitonic hierarchies can be 
generated by (semi) Riemannian metrics Qij on a manifold V of dimension 
diml^ = n > 2 if the geometrical objects are lifted into the total space 
of the tangent bundle TV, or of a vector bundle S = {M,7i,E), dim£^ = 
m > n, by a moving parallel frame formulation of geometric curve flows 
when V has constant matrix curvature as defined canonically with respect to 
certain preferred frames. The second purpose is to elaborate applications in 
geometric mechanics, in particular, that the dynamics defined by any regular 
Lagrangian can be encoded in terms of bi-Hamiltonian structures and related 
solitonic hierarchies. It will be emphasized that, in a similar way, any solution 
of the Einstein equations given by a generically off-diagonal metric can be 
mapped into solitonic equations. 

The paper is organized as follows: 

In section 2 we outline the geometry of vector bundles equipped with a 
nonlinear connection. We emphasize the possibility to define fundamental 
geometric objects induced by a (semi) Riemannian metric on the base space 
when the Riemannian curvature tensor has constant coefficients with respect 
to a preferred nonholonomic framing. 

In section 3 we consider curve flows on nonholonomic vector bundles. 
We sketch an approach to classification of such spaces defined by conven- 
tional horizontal and vertical symmetric (semi) Riemannian subspaces and 
equipped with nonholonomic distributions defined by a nonlinear connection 
structure. Bi-Hamiltonian operators are then derived for a canonical distin- 
guished connection, adapted to the nonlinear connection structure, for which 
the distinguished curvature coefficients are constant. 

Section 4 is devoted to the formalism of distinguished bi-Hamiltonian op- 
erators and vector soliton equations for arbitrary (semi) Riemannian spaces 
admitting nonholonomic deformations to symmetric Riemannian spaces. We 
define the basic equations for nonholonomic curve fiows. Then we consider 
the properties of cosymplectic and symplectic operators adapted to the non- 
linear connection structure. Finally, we construct solitonic hierarchies of 
bi-Hamiltonian anholonomic curve fiows. 

Section 5 contains some applications of the formalism in modern mechan- 
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ics and gravity. 

We conclude with some further remarks in section 6. The Appendix con- 
tains necessary definitions and formulas from the geometry of nonholonomic 
manifolds. 



2 Nonholonomic Structures on Manifolds 

In this section, we prove that for any (semi) Riemannian metric gij on a man- 
ifold V it is possible to define lifts to the tangent bundle TV provided with 
canonical nonlinear connection (in brief, N-connection), Sasaki type metric 
and canonical linear connection structure. The geometric constructions will 
be elaborated in general form for vector bundles. 



2.1 N— connections induced by Riemannian metrics 

Let £ = {E, 71, F, M) be a (smooth) vector bundle of over base manifold M, 
with dimensions dim M = n and dim E = (n+m), for n > 2, and with m > n 
being the dimension of typical fiber F. Here ir : E ^ M defines a surjective 
submersion. At any point u & E, the total space E splits into "horizontal" , 
Mu, and "vertical", Fu, subspaces. We denote the local coordinates in the 
form u = {x,y), or = {x^,y"'), with horizontal indices i,j,k, etc. = 
1, 2, . . . , and vertical indices a, 6, c, etc. = n + l,?7. + 2, ...,n + m|ll The 
summation rule on repeated "upper" and "lower" indices will be applied. 

Let the base manifold M be equipped with a (semi) Riemannian metric, 
namely a second rank tensor of fixed signature|§ hg = g^.{x)dx'^ ® dx^ . It 

is possible to introduce a vertical metric structure vg = g^^{x)dy"' dy^ 

by completing the matrix g..{^) diagonally with ±1 till any nondegenerate 

— u 

second rank tensor g^^{x) li m > n. This defines a metric structure g = 
[hg,vg] (we shall also use the notation g^^ = [g.^,g^i^]) on £. We can deform 
the metric structure, g_^^ — > gaf3 = [gij, gab], by considering a frame (vielbein) 
transformation, 



9at3ix, y) = e^{x, y) e/(x, y)g^{x), (1) 

with coefficients g^^ix) = gapix). The coefficients e^{x,y) will be defined 
below (see formula (JTSll ) from the condition of generating curvature tensors 
with constant coefficients with respect to certain preferred frames. 



^In the particular case when we have a tangent bundle E=TM, both type of indices 
run over the same values since n = m but it is convenient to distinguish the horizontal 
and vertical ones by using different groups of Latin indices. 

^In physics literature, one uses the term (pseudo) Riemannian/Euclidean space 
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For any gab obtained from gap, we consider a generating function 

^x,y) = gah{.x,y)y"'y^ 

inducing a vertical metric 



~9ab = nT^Z;^b (2) 



1 d'^C 

2 dy"-dy 

which is "weakly" regular if det \gab\ 7^ 0. [fl 

By straightforward calculations we can prove the following: 



Theorem 2.1 The Euler-Lagrange equations on TM, 

d fdL\ dL _^ 
dr \dy'J dx' ~ ' 



for the Lagrangian L = y\C\, where = -Jr for a path curve x*(r) on 
M, depending on parameter t, are equivalent to the "nonlinear" geodesic 
equations 

d'^x' , I. dx^ 



defining path curves of a canonical semispray S = y^^ — 2G*(x, y)-^, where 



with g'^ being inverse to the vertical metric 

This theorem has an important geometric mechanical interpretation. 

Corollary 2.1 For any (semi) Riemannian metric g.j{x) on M, we can 
associate canonically an effective regular Lagrange mechanics on TM with the 
Euler-Lagrange equations transformed into nonlinear (semispray) geodesic 
equations. 

The differential of map vr : M is defined by fiber preserving mor- 

phisms of the tangent bundles TE and TM. Its kernel is just the vertical 
subspace vE with a related inclusion mapping i : vE TE. 



■^Similar values, for BcT = 5,^, where (5a— is the Kronecker symbol, were introduced for 
the so-called generalized Lagrange spaces when C was called the "absolute energy" [2] . 
^See Refs. [H |3] for details of a similar proof; here we note that in our case, in general, 
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Definition 2.1 A nonlinear connection (N- connection) on a vector bun- 
dle £ is defined by the splitting on the left of an exact sequence 



^ vE ^ TE ^ TE/vE 0, 

i.e. by a morphism of submanifolds N : TE vE such that Noi is the 
identity map in vE. 

Equivalently, a N-connection is defined by a Whitney sum of conventional 
liorizontal (li) subspace, (hE), and vertical (v) subspace, (vE), 

TE = hE® vE. (3) 

This sum defines a nonholonomic (alternatively, anholonomic, or non-integ- 
rable) distribution of horizontal (h) and vertical (v) subspaces on TE. Lo- 
cally, a N-connection is determined by its coefficients A^f (m), 

N = N;'{u)dx'® ^ 



The well known class of linear connections consists of a particular subclass 



with the coefficients being linear in y'^, i.e., N°-{u) = TlAx)y'^. 



Remark 2.1 A manifold (or a bundle space) is called nonholonomic if it is 
provided with a nonholonomic distribution (see historical details and sum- 
mary of results in l2^). In a particular case, when the nonholonomic distri- 
bution is of type such spaces are called N-anholonomic f2t 



Any N-connection N = {N°-{u)} may be characterized by a N-adapted 
frame (vielbein) structure e^ = (ei,ea), where 

and the dual frame (coframe) structure e'^ = (e*,e°), where 

e' = dx' and e'^ = dy" + N^{u)dx\ (5) 

We remark that = (ej, Ca) and e'^ = (e*, e") are, respectively, the former 
"N-elongated" partial derivatives 6^, = = {Si, da) and N-elongated 

differentials S'^ = 5u^ = {d'^,S°') which emphasize that operators (jlj) and 
([5]) define, correspondingly, certain "N-elongated" partial derivatives and 
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differentials which are more convenient for tensor and integral calculations 
on such nonholonomic manifoldsj^l 

For any N-connection, we can introduce its N-connection curvature 

n = In^. d' A d' ® da, 

2 ■' 

with the coefficients defined as the Nijenhuis tensor, 

V i\ 3 t t 3 Q^j Q^i * Qyh 3 Qyh ^ > 

The vielbeins ^ satisfy the nonholonomy (equivalently, anholonomy) 
relations 

[e^,, ep] = e^efs - efse^ = W^^^^ (7) 

with (antisymmetric) nontrivial anholonomy coefficients = daN^ and 

These geometric objects can be defined in a form adapted to a N-connec- 
tion structure via decompositions being invariant under parallel transport 
preserving the splitting In this case we call them "distinguished" (by 
the N-connection structure), i.e. d-objects. For instance, a vector field 
X G TV is expressed 

X = {hX, vX), or X = X"e« = X^e^ + X'^e„, 

where hX = X*ej and vX = X'^Ca state, respectively, the N-adapted hori- 
zontal (h) and vertical (v) components of the vector (which following Refs. 
[2l[3j is called a distinguished vector, in brief, d-vector). In a similar fashion, 
the geometric objects on V, for instance, tensors, spinors, connections, etc. 
can be defined and called respectively d-tensors, d-spinors, d-connections if 
they are adapted to the N-connection splitting 

Theorem 2.2 Any (semi) Riemannian metric g-j{x) on M induces a canon- 
ical N-connection on TM. 

Proof. We sketch a proof by defining the coefficients of N-connection 

iV",(.,!/)=g- (8) 



^We shall use "boldface" symbols whenever necessary as emphasis for any space and/or 
geometrical objects equipped with/adapted to a N-connection structure, or for the coef- 
ficients computed with respect to N-adapted frames. 
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where 

fim = ^9'''{dmgih + digmh-dhgim), dh = d/dx^, 

with gah and c/ij defined respectively by formulas ^ and (j2]). □ 

The N-adapted operators (jl]) and ([5]) defined by the N-connection coef- 
ficients dH]) are denoted respectively e,^ = (ei,ea) and e'^ = (e*,e"). 

2.2 Canonical linear connection and metric structures 

The constructions will be performed on a vector bundle E provided with N- 
connection structure. We shall emphasize the special properties of a tangent 
bundle {TM,Tr,M) when the linear connection and metric are induced by a 
(semi) Riemannian metric on M. 

Definition 2.2 A distinguished connection (i.e., d-connection) D = {hD, 
vJD) is a linear connection preserving under parallel transport the nonholo- 
nomic decomposition 

The N-adapted components F"^^ of a d-connection Dq = (e^jD) are 
defined by equations 

D,e^ = r^^^e^, or T\p (u) = (D,e^)Je^. (10) 

The N-adapted splitting into h- and v-covariant derivatives is stated by 

hB = {D, = and vB = {D, = {q„C^,)}, 

where, by definition, L}^ = (Dfcej)Je\ L^^ = (Dfceb)Je", Cj^ = (D^ejOJeS 
= (D,eb)Je^ The components = {L]„ Lt„q^,C^^) completely 

define a d-connection D on E. 

The simplest way to perform N-adapted computations is to use differen- 
tial forms. For instance, starting with the d-connection 1-form, 

F°^ = F'^^^e^ (11) 

with the coefficients defined with respect to N-elongated frames (I5]) and (jll), 
the torsion of a d-connection, 

= De° = rfe" + F"^ A e^, (12) 
8 



is characterized by (N-adapted) d-torsion components, 



rpl TZ TZ rpl rpl j^z rpa Qa 

jk ~ jk kj^ ^ ja ~ ^ aj ~ ^ ja^ ^ ji ~ ji^ 

dm 

rpa rpa ~_ I ja rpa r^a r^a ^^'X\ 

^ hi — -'- ib — ^ bi^ be — ^ be ^ cb- \^'^) 

For d-connection structures on TM, we have to identify indices in the form 
i ^ a, j ^ b, ... and the components of N- and d-connections, for instance, 
^ and L*^-^ L\i., ^ C^^^ ^ C-ji^. 

Definition 2.3 A distinguished metric (i.e., d-metric) on a vector bundle 
El is a non- degenerate second rank metric tensor g =g(B]sih, equivalently 

g = gij{x, y) e ® + hab{x, y) e" ® e^ (14) 

adapted to the N-connection decomposition 

From the class of arbitrary d-connections D on V, one distinguishes those 
which are metric compatible (metrical) satisfying the condition 

Dg = (15) 

including all h- and v-projections Djgki = 0, Dagn = 0, Djhab = 0, Dahhc = 

0. For d-metric structures on V ~TM, with gij = hab, the condition of 
vanishing "nonmetricity" (ITSl) transforms into 

HBig) =0 and vBih) =0, (16) 

1. e. Djgki = and Dagn = 0. 

For any metric structure g on a manifold, there is the unique metric com- 
patible and torsionless Levi Civita connection V for which Vg = 0. This con- 
nection is not a d-connection because it does not preserve under parallelism 
the N-connection splitting One has to consider less constrained cases, 
admitting nonzero torsion coefficients, when a d-connection is constructed 
canonically for a d-metric structure. A simple minimal metric compatible 
extension of V is that of canonical d-connection with T^^, = and T°^^ = 
but T^ja, T'^ji and T\ are not zero, see (fT^ . The coefficient formulas for such 
connections are given in Appendix, see (lU^ and related discussion. 

Lemma 2.1 Any (semi) Riemannian metric g.^{x) on a manifold M induces 
a canonical d-metric structure on TM, 

g = 9ij{x, y) (g) + gzj{x, y) e' ^ & , (17) 
where e* are elongated as in but with N^j given by 
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Proof. This construction is trivial by lifting to the so-called Sasaki 
metric [2Zj but in our case using the coefficients (jij I^^.D 

Proposition 2.1 There exist canonical d- connections on TM induced by a 
(semi) Riemannian metric g--{x) on M . 

Proof. We can construct an example in explicit form by introducing gij 
and Qab in formulas flTUj) . see Appendix, in order to compute the coefficients 

The above Lemma and Proposition establish the following result. 

Theorem 2.3 Any (semi) Riemannian metric g--{x) on M induces a non- 
holonomic (semi) Riemannian structure on TM. 

We note that the induced Riemannian structure is nonholonomic because 
on TM there is a non-integrable distribution ([H]) defining N^-. The cor- 
responding curvature curvature tensor R'^pr^^. = {R\jk: P^jkai ^\cd} com- 
puted by substituting cjij, N^j and e^ into formulas fl75l) . from Appendix, for 

= (L^jky^bc)- Here one should be noted that the constructions on TM 
depend on arbitrary vielbein coefficients e^{x,y) in ([T]). We can restrict 
such coefficients in order to generate various particular classes of (semi) Rie- 
mannian geometries on TM, for instance, in order to generate symmetric 
Riemannian spaces with constant curvature, see Refs. [HI [2H1 ES] • 

Corollary 2.2 There are lifts of a (semi) Riemannian metric g.j{x) on M, 
dimM = n, generating a Riemannian structure on TM with the curva- 
ture coefficients of the canonical d-connection coinciding (with respect to N- 
adapted bases) with those for a Riemannian space of constant curvature of 
dimension n + n. 

Proof. For a given metric g..{x) on M, we chose such coefficients e^{x, y) 
= {ej^{x,y)} in ([I]) that 

9ab{x,y) = e^-{x,y) ef{x,y)gjyx) 

produces a vertical metric ([2]) of type 

= 2d^ = 2 drdyf ^^^^^ = ^'^^ 

where gab is the metric of a symmetric Riemannian space (of constant cur- 
vature). Considering a prescribed gab, we have to integrate two times on 
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in order to find any solution for e^- defining a frame structure in tlie vertical 
subspace. Tfie next step is to construct tlie d-metric c/a/s = [ gij, cjab] of 
type (fTTl) . in our case, witli respect to a nonfiolonomic base elongated by 

N\, generated by g..{x) and get = gab, like in and IQ. This defines a 
constant curvature Riemannian space of dimension n + n. The coefficients 
of the canonical d-connection, which in this case coincide with those for the 
Levi Civita connection, and the coefficients of the Riemannian curvature can 
be computed respectively by putting g^f = gab in formulas (1701) and fl73|) . see 
Appendix. Finally, we note that the induced symmetric Riemannian space 
contains additional geometric structures like the N-connection and anholon- 
omy coefficients W^p, see (I7]).n 

There are various possibilities to generate on TM nonholonomic Rieman- 
nian structures from a given metric g^-{x) on M. They result in different 
geometrical and physical models. In this work, we emphasize the possibil- 
ity of generating spaces with constant curvature because for such symmetric 
spaces one can derive a bi-Hamiltonian hierarchy of curve flows and associ- 
ated solitonic equations. 

Example 2.1 The simplest example where a Riemannian structure with con- 
stant matrix curvature coefficients is generated on TM comes from a d- 
metric induced by gij = 6ij, i.e. 

SlE] = Sijc' ® + 6ij ® (19) 

with e* given by N^j as defined by g^.{x) on M. 

It should be noted that the metric ( |T9|) is generically off-diagonal with 
respect to a coordinate basis because, in general, the anholonomy coefficients 
from ([7]) are not zero. This way, we model on TM a nonholonomic Euclidean 
space with vanishing curvature coefficients of the canonical d-connection (it 
can be verified by putting respectively the constant coefficients of metric ( fT9l) 
into formulas (|70l) and ( |75l) ). We note that the conditions of Theorem 12. II are 
not satisfied by the d-metric (|T9l) (the coefficients gij = Sij are not defined 
as in ([2])), so we can not associate a geometrical mechanics model for such 
constructions. 

There is an important generalization: 

Example 2.2 We can consider C as a hypersurface in TM for which the 
matrix d'^C/dy"'dy^ (i.e. the Hessian, following the analogy with Lagrange 
mechanics and field theory) is constant and nondegenerate. This states that 
cjij = const, which produces vanishing curvature coefficients for the canonical 
d-connection induced by g^.{x) on M . 
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Finally, we note that a number of geometric ideas and methods applied 
in this section were considered in the approaches to the geometry of non- 
holonomic spaces and generalized Finsler-Lagrange geometry elaborated by 
the schools of G. Vranceanu and R. Miron and by A. Bejancu in Romania 
[301 [3T| El El ESI El] , see also Section [5l In these approaches it is possible to 
construct geometric models with metric compatible linear connections, which 
is important for developments connected with modern (non) commutative 
gravity and string theory [25l EE]. For Finsler spaces with nontrivial non- 
metricity, for instance, those defined by the the Berwald and Chern connec- 
tions [32], the resulting physical theories with local anisotropy fall outside of 
the class of standard models. 

3 Curve Flows and Anholonomic Constraints 

We now formulate the geometry of curve flows adapted to the nonlinear 
connection structure. 

3.1 N— adapted curve flows 

Let us consider a vector bundle £ = {E,!!, F, M), dimE = n + m (in the 
particular case, E = TM, we have m = n) equipped with a d-metric g = 
[g, h] (|T^ and N-connection N°- structures. A non-stretching curve 7(r, /) 
on V, where r is a parameter and / is the arclength of the curve on V, is 
defined by an evolution d-vector Y = and tangent d-vector X = 7; such 
that g(X, X) = 1. Such curves 'y{T,l) sweep out a two-dimensional surface 
in T^(r,z)V C TV. 

We shall work with N-adapted bases (jl]) and and the connection 1- 
form r"g = r^g^e''' with the coefficients r°g^ for the canonical d-connection 
operator D ( l69l) (see Appendix), acting in the form 

Dxe, = (Xjrj)e^ and Dve, = (YjrJ)e^, (20) 

where "J " denotes the interior product and the indices are lowered and raised 
respectively by the d-metric ga/s = [gij,hab] and its inverse g°'^ = [g''^,h°'^]. 
We note that Dx = X°D„ is the covariant derivation operator along curve 
7(r, /). It is convenient to fix the N-adapted frame to be parallel to curve 
7(Z) adapted in the form 

= hX, for z = 1, and e^, where hg{hX,e^) =0, (21) 
e"'+^ = vX, for a = n + 1, and e", where vg{vX, e") =0, 
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for i = 2,3, ... ,n and a = n + 2,n + 3, . . . ,n + m. For such frames, the 
covariant derivative of each "normal" d-vector e" is parallel to the d-vectors 
adapted to 7(r, /), 



(22) 



Dxe' = -p'{u) X and Bh^hX = p'{u) e., 
Dxe^ = -p^{u) X and D^x^^X = p^{u) e^, 

in terms of some coefficient functions and p°'{u). The formulas 

and (|22|) are distinguished into h- and v-components for X =hX + vX and 
D = [hJ^, vT>) for D = {r^,^}, hT> = {L]^, L^J and t;D = {Cj,, C^]. Along 
7(/), we can pull back differential forms in a parallel N-adapted form. For 
instance, r";J = Xjr"^. 

An algebraic characterization of parallel frames can be obtained if we 
perform a frame transformation preserving the decomposition ^ to an or- 
thonormal basis e.^', 

e^^A^\u)e^,, (23) 
called an orthonormal d-basis, where the coefficients of the d-metric ffT^ are 



transformed into the Euclidean ones ga'/j' 
obtain two skew matrices 



6. 



a'l3'- 



In distinguished form, we 



= hX\r^ = 2 e^^ p>'^ and V^J = vX\r- = 2 e^x p' 



Wb' 



where 



g{hX,e'') = [1, 0, . . . , 01 and e^x ^ h{vX,e''') = [1,0,... 



m—l 



and 





-Pi' 



On 



and r,x./' 





-pa 



P 
Or, 



with 0[h] and 0[y] being respectively (n — 1) x (n — 1) and (m — 1) x (m — 1) 
matrices. The above presented row-matrices and skew-matrices show that 
locally a N-anholonomic manifold V of dimension n + m, with respect to 
distinguished orthonormal frames, are characterized algebraically by pairs 
of unit vectors in and preserved respectively by the SO{n — 1) and 
SO{m — 1) rotation subgroups of the local N-adapted frame structure group 
SO{n)® SO{m). The connection matrices F^^i' and F^^a' ^' belong to the 
orthogonal complements of the corresponding Lie subalgebras and algebras, 
so{n — 1) C so{n) and 5o(m — 1) C so(m). The torsion ( JT2l) and curvature 
(ITT]) (see Appendix) tensors in orthonormal component form with respect to 
(|2T|) can be mapped into a distinguished orthonormal dual frame ( 123|) . 



= Dxe^ - DveS: 



t'Y-'^ X/3' 



(24) 
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and 



(X, Y) - DyFx^" - DxFy^" + Ty^T Fx^" - F^^T Fy^? , (25) 

where By = g(Y,e"') and Fy^"' = YJF^' = g(e"', Dve/j') define re- 
spectively the N-adapted orthonormal frame row-matrix and the canoni- 
cal d-connection skew-matrix in the flow directions, where TZ^, (X, Y) = 
g(e"', [Dx, DY]e^/) is the curvature matrix. Both torsion and curvature 
components can be distinguished into h- and v-components like ( |T3l) and 
(172|) . by considering N-adapted decompositions of type 

g = [g, h],e,3, = (e^v, Cb'), e°' = (e^', e"'), X =/iX + i;X, D = (/iD, vB). 

Finally, we note that the matrices for torsion (!24l) and curvature (|25ll can be 
computed for any metric compatible linear connection like the Levi Civita 
and the canonical d-connection. For our purposes, we are interested to define 
a frame in which the curvature tensor has constant coefficients and the torsion 
tensor vanishes. 

3.2 Anholonomic bundles with constant matrix curva- 
ture 

For vanishing N-connection curvature and torsion, we get a holonomic Rie- 
mannian manifold with constant curvature, wherein the equations fl2^ and 
( !25l) directly encode a bi-Hamiltonian structure, see details in Refs. [T3l [T^ . 
A larger class of Riemannian manifolds for which the frame curvature ma- 
trix is constant consists of the symmetric spaces M = G/H for compact 
semisimple Lie groups G D H (where H is required to be invariant under 
an involutive automorphism of G). A complete classification and summary 
of main results for such spaces is given in Refs. fl9[ [28] . Constancy of the 
frame curvature matrix is a consequence of the fact that the Riemannian 
curvature and the metric tensors on the curved manifold M = G/H are co- 
variantly constant and G-invariant. In Ref. [Uj, a bi-Hamiltonian structure 
was shown to be encoded in the frame equations analogous to fl2^ and fl23|) 
for the symmetric spaces M = G/SO{n) with H = SO{n) D 0{n — 1). 
All such spaces are exhausted by G = SO{n + 1), SU{n). The derivation of 
the bi-Hamiltonian structure exploited the canonical soldering of the spaces 
G/SO{n) onto Klein geometries [29j. A similar derivation [33j holds for the 
Lie groups G = SO{n + 1), SU{n) themselves when they are viewed as sym- 
metric spaces in the standard manner [2H]- A broad generalization of this 
derivation to arbitrary Riemannian symmetric spaces, including arbitrary 
compact semisimple Lie groups, has been obtained in Ref. [31] . 
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3.2.1 Symmetric nonholonomic tangent bundles 

We suppose that the base manifold is a constant-curvature symmetric space 
of the form M = hG / SO{n) with the isotropy subgroup hH = SO{n) D 0(n) 
and that the typical fiber space is likewise of the constant-curvature symmet- 
ric form F = vG / SO{m) with the isotropy subgroup vH = SO{m) D 0{m). 
This means (according to Cartan's classification of symmetric spaces |19j ) 
that we have hG = S0{n + 1) and vG = S0{m + 1), which is general enough 
for a study of real holonomic and nonholonomic manifolds and geometric 
mechanics models. 

Our aim is to solder in a canonical way (like in the N-connection geome- 
try) the horizontal and vertical symmetric Riemannian spaces of dimension 
n and m with a (total) symmetric Riemannian space V of dimension n + m, 
when V = G/ SO{n+m) with the isotropy group H = SO{n+m) D 0{n+m) 
and G = SO{n + m + 1). First, we note that the just mentioned horizon- 
tal, vertical and total symmetric Riemannian spaces can be identified with 
respective Klein geometries. For instance, the metric tensor hg = {cjij} on 
M is defined by the Cartan-Killing inner product < ■, ■ >h on T^hG ^ hg 
restricted to the Lie algebra quotient spaces hp =hg/h[), with T^hH ^ hi), 
where hg =h\) © hp is stated such that there is an involutive automorphism 
of hG leaving hH fixed, i.e. [/;,f),/;,p] C hp and [/ip,/ip] C h[). In a similar 
form, we can define the inner products and Lie algebras: vg = {hat} is given 
by restriction of < ■,■ >t, on TyvG ~ vq to vp =vq/v[) with TyvH ~ v\), 
VQ =f f) © f p, where 

[v\),vp\ C ^;p, [vp,vp\ C v\)- (26) 

likewise g = {gajs} is given by restriction of < ■, ■ >g on T(^x,y)G ~ to 
p = with T(x,y)H ~ P), g = [) © p where 

We parametrize the metric structure with constant coefficients on V = 
G/SO{n + m) in the form 

g = gapdu"' © du^, 
where are local coordinates and 



9al3 



h^ + NlN^^ha, N^hae 
Nfhbe hah 



(27) 



^It is necessary to consider hG = SU{n) and vG — SU{m) for geometric models 
involving spinor and gauge fields 
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where trivial, constant, N-connection coefficients iVJ = h Qj^ are defined 

in terms of any given coefficients h"^^ and cjjh-, i.e. from the inverse metrics 
coefficients defined respectively by the metric on hG = SO{n + 1) and by 
the metric coefficients out of {n x n)- and (m x m)- terms of the metric cjap. 
As a result, we define an equivalent d-metric structure of type (fHIl 

g = Qi, e ® + Kb ® e^ (28) 

defining a trivial (n + m)-splitting g =g(B-^h because all nonholonomy co- 
efficients WCfl and N-connection curvature coefficients Cl% are zero. In a 
more general form, we can consider any covariant coordinate transforms of 
( l28l) preserving the (n + m)-splitting resulting in any W^p = (C!) and 
= (E]). We denote such trivial N-anholonomic Riemannian spaces as 
il =[hG = SO{n + 1), vG = SO{m + l), iVf]. It can be considered that 
such trivially N-anholonomic group spaces possess a Lie d-algebra symmetry 
so^(n + m) =F so(n) © so(m). 

The simplest generalization for a vector bundle E is to consider non- 
holonomic distributions on = G/SO{n + m) defined locally by arbitrary 
N-connection coefficients N°'{x,y) with nonvanishing W^/^ and but with 
constant d-metric coefficients when 

g = e' ® + hab ® e^ (29) 
= dx\ e" = dy" + N^{x, y)dx\ 

This metric is very similar to (|T9l) but with the coefficients cjij and hab in- 
duced by the corresponding Lie d-algebra structure S0f:j{n + m). Such spaces 
transform into N-anholonomic Riemann-Cartan manifolds Vn = [hG = 
SO{n + 1), vG = SO{m + 1), N^] with nontrivial N-connection curvature 
and induced d-torsion coefficients of the canonical d-connection (see for- 
mulas (fT3ll computed for constant d-metric coefficients and the canonical 
d-connection coefficients in (|69l) ). One has zero curvature for the canoni- 
cal d-connection (in general, such spaces are curved ones with generically 
off-diagonal metric fl^IH) and nonzero curvature tensor for the Levi Civita 
connection) This allows us to classify the N-anholonomic manifolds (and 
vector bundles) as having the same group and algebraic properties as pairs 
of symmetric Riemannian spaces of dimension n and m but nonholonomi- 
cally soldered to the symmetric Riemannian space of dimension n + m. With 



^With constant values for the d-metric coefficients we get zero coefficients for the 
canonical d-connection which then yields vanishing curvature coefficients (|72p. 
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respect to N-adapted orthonormal bases with distinguished h- and v- 
subspaces, we obtain the same inner products and group and Lie algebra 
spaces as in ( |26l) . 

The classification of N-anholonomic vector bundles is almost similar to 
that for symmetric Riemannian spaces if we consider that n = m and try 
to model tangent bundles of such spaces, provided with N-connection struc- 
ture. For instance, we can take a (semi) Riemannian structure with the 
N-connection induced by a absolute energy structure like in ([8]) and with the 
canonical d-connection structure (l69ll . for g^f = gab ([18]) • A straightforward 
computation of the canonical d-connection coefficientqj and of d-curvatures 
for °gij and °A^*j proves that the nonholonomic Riemannian manifold 
(M = SO{n + 1)/S0{n), °C) possess constant both zero canonical d-con- 
nection curvature and torsion but with induced nontrivial N-connection cur- 
vature Such spaces, being tangent to symmetric Riemannian spaces, 
are classified similarly to the Riemannian ones with constant matrix curva- 
ture, see ( !26|) for n = m, but possessing a nonholonomic structure induced 
by generating function °C 

3.2.2 anholonomic Klein spaces 

The bi-Hamiltonian and solitonic constructions in Refs. [SU [T71 dU are 
based on soldering Riemannian symmetric-space geometries onto Klein ge- 
ometries [29]. For the N-anholonomic spaces of dimension n + n, with con- 
stant d-curvatures, similar constructions will hold but we have to adapt them 
to the N-connection structure. 

There are two Hamiltonian variables given by the principal normals 
and ^u, respectively, in the horizontal and vertical subspaces, defined by the 
canonical d-connection D = (/iD,uD), see formulas ( 12T1) and ( l22l) . 

= D;,x/iX = iy% and V = D^x^^X = u^e^. 

This normal d-vector v = ( ^u, ^u), with components of type i/" = (z/*, 
u"-) = (z/^, z/*, z/*^"^^, z/"), is defined with respect to the tangent direction of 
curve 7. There is also the principal normal d-vector w = ( ^zu, ^w) with 
components of type ro" = {w\ zu"') = {zu^ , zu\ zu"'^^ , w"') defined with 
respect to the flow direction, with 

'^zu = D/iY^X =ro*ej, "zo = Dt,YfX = tzr^ea, 
**0n tangent bundles, such d-connections can be defined to be torsionless 
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representing a Hamiltonian d-covector field. We can consider that the normal 
part of the flow d-vector 

represents a Hamiltonian d-vector field. For such configurations, we can 
consider parallel N-adapted frames e^/ = (ej/,ea') when the h-variables z/*', 
zu^',h^' are respectively encoded in the top row of the horizontal canonical 
d-connection matrices F^^i' ^hYi' matrix (gy)^ =f 

Gy — g\\ Gx where = g(hY,h'X.) is the tangential h-part of the flow d- 
vector. A similar encoding holds for v-variables u"-' , w""' , h"-' in the top row 
of the vertical canonical d-connection matrices T^^a' ^ ^vYa' ^ 
the row matrix (eY)_|_ == Gy — /iy where = h{vY,vX.) is the tangential 
v-part of the flow d-vector. In an abbreviated notation, we shall write 
v" and Ct7° where the primed Greek indices .•• will denote both N- 

adapted and then orthonormalized components of geometric objects (like d- 
vectors, d-covectors, d-tensors, d-groups, d-algebras, d-matrices) admitting 
further decompositions into h- and v-components defined as non-integrable 
distributions of such objects. 

With respect to N-adapted orthonormal frames, the geometry of N- 
anholonomic manifolds is defined algebraically, on their tangent bundles, by 
pairs of horizontal and vertical Klein spaces (see [29j for a summary of Klein 
geometry, and for bi-Hamiltonian soliton constructions, see [Ml [Hj). The 
N-connection structure induces a N-anholonomic Klein geometry stated by 
two left-invariant hg- and f0-valued Maurer-Cartan forms on the Lie d- 
group G = [hG, vG) that are identified with the zero-curvature canonical 
d-connection 1-form = { *^r"^,}, via 

G-pa' _ G-pa' W _ hG j i' i vG/^i' ^k' 

For the case of a trivial N-connection structure in vector bundles with the 
base and typical fiber spaces being symmetric Riemannian spaces, we identify 
^^U-iy and ^^C^jif^, with the coefficients of the Cartan connections ^^L 
and ^^C, respectively, for the hG and vG, both with vanishing curvatures. 
We have 

rf^r+i[^r, ^r] = o 

and, respectively, for h- and v-components, d + ^['^'^L, '^'^L] = 
and d + |[ "'^C, ""^C] = 0, where d denotes the total derivative on the 
d-group manifold G = hG(BvG or its restrictions on hG or vG. We can con- 
sider that defines the so-called Cartan d-connection for non-integrable 
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N-connection structures (see details and supersymmetric/ noncommutative 
developments in [2Sl[2n])- 

The Cartan d-connection determines a N-anholonomic Riemannian geo- 
metric structure on the nonholonomic bundle 

E=[hG = SO{n + 1), vG = SO{m + 1), N^], 

derived through the following decomposition of the Lie d-algebra q = hg(Bvg. 
For the horizontal splitting, we have hg = so(n) © hp, with [hp, hp] C so(n) 
and, for the vertical splitting, we have [so{n),hp] C hp; for the vertical 
splitting, vg = so{m)(Bvp, with [fp,fp] C so(m) and [so{m),vp] C vp. When 
n = m, any canonical d-objects (N-connection, d-metric, d-connection, etc.) 
derived from (|29l) using the Cartan d-connection agree with ones determined 
by the canonical d-connection (!70l) on a tangent bundle. 

It is useful to consider a quotient space with distinguished structure group 
Vn = G/ SO{n)Q) SO{m) regarding G as a principal {SO{n) ® SO{m))- 
bundle over E, which is a N-anholonomic bundle. In this case, we can always 
fix a local section of this bundle and pull-back to give a {hg © vg)- 
valued 1-form in a point m G E. Any change of local sections define 
SO{n)® SO{m) gauge transformations of the canonical d-connection ^V, all 
preserving the nonholonomic decomposition 

There are involutive automorphisms ha = ±1 and va = ±1, respectively, 
of hg and vg, defined that 50(77.) (or so(m)) is the eigenspace ha = +1 (or 
va = +1) and hp (or vp) is the eigenspace ha = —1 (or va = —1). By means 
of an N-adapted decomposition relative to these eigenspaces, the symmetric 
part of the canonical d-connection 

r=l(«r+a(«r)), 

defines a (so(n) © 5o(m))-valued d-connection 1-form, with respective h- 
and v-splitting L =1 {^^L+ha {^^L)) and C =i(^9C + hai^^C)). Likewise 
the antisymmetric part 

e=l(9r-a(T)), 

defines a {hp © w p)-valued N-adapted coframe for the Cartan-Killing inner 
product < ■, ■ >p on T„G ~ hg Q) vg restricted to T^Vn — p, with respec- 
tive h- and v-splitting he =1 {''^L-ha C^sL)) and ve =^("^0 - hai^"^C)). 
This inner product, distinguished into h- and v-components, provides a 
d-metric structure of type g = [g,h] ( HM . where g =< he®he >hp and 
h =< veiS)ve >yp on Vn = G/S'0(n)© SO{m). 
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We generate a G-invariant d-derivative D whose restriction to the tan- 
gent space TVn for any N-anhofonomic curve flow 7(r, /) in Vn = G/ SO{n) 
© SO{rn) is defined via 

Dxe = [e,7jr] and Dve = [e, 7,jr] , (30) 

which have h- and v-decompositions. The derivatives Dx and Dy are equiv- 
alent to those considered in ( l20l) and obey the Cartan structure equations 
( 12^ and ( l25i) . For canonical d-connections, a large class of N-anholonomic 
spaces of dimension n = m, the d-torsions are zero and the d-curvatures 
have constant coefficients. 

Let e° = (e* , e" ) be a N-adapted orthonormal coframe being identi- 
fied with the {hp © f p)-valued coframe e in a fixed orthonormal basis for 
p =hp © vp Ghg © vg. Considering the kernel/ cokernel of Lie algebra mul- 
tiplications in the h- and v-subspaces, respectively, [e/^x, ■]hQ [e^x, 
we can decompose the coframes into parallel and perpendicular parts with 
respect to ex- We write 

e = (e,^ = hec + vec, e^x = hec± + f ec±), 

for p( = hp Q) t>p)-valued mutually orthogonal d- vectors ec and e^-x, sat- 
isfying the conditions [exjecj^ = and [ex,ecx]g 7^ 0; such conditions 
can be stated in h- and v-component form, respectively, [he^, hec]j^g = 0, 
[/iex,^ec.x]^g 7^ and [vex,vec]^^ = 0, [vex,vec±]^^ ^ 0. Then we have 
the algebraic decompositions 

T„Vn ^ p =hp ®vp = Q =hQ © vg/5o{n) © so(m) 

and 

P = Pc®Pc^ = {hpc ® vpc) © {hpc± © vpc±) , 

with p|| C pc and pc± C p_L, where [p||,pc] = 0, < pc^,Pc >= 0, but 
[p||, pcx] 7^ (i.e. pc is the centralizer of ex in p =hp © f p chg © vg); in h- 
and v-components, thus hp\\ C hpc and hpc^- C hp±, where [/ipy, /ipc] = 
0, < hpcJ-,hpc >= 0, but [/ip||, /ipcx] 7^ (i.e. hpc is the centralizer 
of e^x in hp <Zhg) and vp\\ C vpc and vpc-L C vp_i, where [wp||,fpc7] = 
0, < fp(7x,t>pc >= 0, but [t>p||, t>pf7x] 7^ (i.e. vpc is the centralizer of 
e^x in ^^P Cfg). Using the canonical d-connection derivative Dx of a d- 
co vector algebraically perpendicular (or parallel) to ex, we get a new d- vector 
which is algebraically parallel (or perpendicular) to ex, i.e. Dxec G p^x (or 
OxGc^ £ Pc)] in h- and v-components such formulas are written Dhxhec G 
hpc± (or Dftx/iecx G hpc) and D^x^^ec G vpc± (or D^x^^ec-x G vpc)- AH 
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such d-algebraic relations can be written in N-anholonomic manifolds and 
canonical d-connection settings, for instance, using certain relations of type 

Dx(e"')c = v";,(e^%x and Dx(e-%x = -v-;,(e'^')c, 

for some antisymmetric d-tensors v°'^' = — v'^'"'. Note we will get a N- 
adapted {SO{n) © S'0(m))-parallel frame defining a generalization of the 
concept of Riemannian parallel frame [13] on N-adapted manifolds whenever 
pc is larger than py. Substituting e" = (e* , ) into the last formulas and 
considering h- and v-components, we define 5'0(n)-parallel and SO{m)- 
parallel frames (for simplicity we omit these formulas when the Greek letter 
indices are split into Latin letter h- and v-indices). 

The final conclusion of this section is that the Cartan structure equations 
on hypersurfaces swept out by nonholonomic curve flows on N-anholonomic 
spaces with constant matrix curvature for the canonical d-connection ge- 
ometrically encode separate 0{n — 1)- and 0{m — l)-invariant (respec- 
tively, horizontal and vertical) bi-Hamiltonian operators, whenever the N- 
connection freedom of the d-group action SO{n)Q)SO{m) on e and F is used 
to fix them to be a N-adapted parallel coframe and its associated canoni- 
cal d-connection 1-form, related to the canonical covariant derivative on 
N-anholonomic manifolds. 

4 Anholonomic bi-Hamiltonian Structures 
and Vector Soliton Equations 

Introducing N-adapted orthonormal bases, for N-anholonomic spaces of di- 
mension n + n, with constant curvatures of the canonical d-connection, we 
will be able to derive bi-Hamiltonian and vector soliton structures similarly 
to the case of symmetric Riemannian spaces (see Ref. [2TI [TTl [M]). 

4.1 Basic equations for N— anholonomic curve flows 

In this section, we shall prove the results for the h-components of certain 
N-anholonomic manifolds with constant d-curvature and then duplicate the 
formulas for the v-components omitting similar details. 

There is an isomorphism between the real space so(n) and the Lie algebra 
of n X n skew-symmetric matrices. This yields an isomorphism between hp 
~ M" and the tangent spaces T^M = so{n + 1)/ so(n) of the Riemannian 
manifold M = SO{n + 1)/ SO{n) as described by the following canonical 
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decomposition 

hQ = so{n + 1) D hp e 





-hp^ 



hp 
hO 



for hO ehi) — so{n) 



with hp = {j5*'} GM"" being the h-component of the d-vector p = 

and hp^ denoting the transposition of the row hp. The Cartan-KiUing inner 

product on hp is given by 



hp-hp 



hp ' 




hp ' 


-hp^ hO 




-hp^ hO 





hp ' 


T 


hp ' 




{ 


_ -hp^ hO _ 




_ -hp^ hO _ 


) 



= -tr 
2 



where tr denotes the trace of the corresponding product of matrices. This 
product identifies canonically hp ~ with its dual hp* ~ R". In a similar 
fashion, we can define the Cartan-KiUing inner product vp-vp where 



VQ = so(m + 1) D e 





7 

—vp 



vp 
vO 



for vO Gwf) = so{m) 



with vp — {p"'} eR"^ being the v-component of the d-vector p = (p* p"')- In 
general, in the tangent bundle of a N-anholonomic manifold, we can consider 
the Cartan-Killing N-adapted inner product p ■ p =hp-hp+vp-vp. 

Following the introduced Cartan-Killing parametrizations, we analyze 
the flow 7(t, /) of a non-stretching curve in Vn = G/ SO{n)Q) SO{m). 
Let us introduce a coframe e e T*Vn ® {hp®vp), which is a N-adapted 
(S'0(n)®5'0(m))-parallcl covcctor frame along 7, and its associated canon- 
ical d-connection 1-form F G T^Vn (E) (5o(n)©5o(m)). Such d-objects are 
respectively parametrized: 



ex = e/ix + e„x, 



for 



and 



e/ix = 7/ixJ/ie 





-(i,"o)^ 



(1,0) 
hO 



e,;x = 7i;xjt'e = 
where we write (1, "O^) e M"*, "O^ e 



(1,0) 





.-(1,V 

and (1, ) e 

[r/ix, r^x] , 



e 



pm— 1. 
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with 



where 



(0, 
-(0,"0r L 



— V 



V 
hO 



e so{n), Ij e 



on— 1 



e 50(n + 1), 



/lO eso{n - 1), 



and also 



where 



^T.xJC 



V 



(0, ) 
-(0,0"f c 



e 5o(m + 1), 



eso(m), v e , vOeso(m-l). 



The above parametrizations involve no loss of generality, because they can be 
achieved by the available gauge freedom using SO{n) and SO{m) rotations 
on the N-adapted coframe and canonical d-connection 1-form, distinguished 

in h- and v-components. 

We introduce decompositions of horizontal SO{n + 1)/ SO{n) matrices 



hp3 



/ip 



+ 





/ip||,"0 


-(0,/ip±)' 



(/ipil,"0 
(0,/ip±) 



into tangential and normal parts relative to e/ix via corresponding decom- 
positions of h-vectors /ip = (/ip||,/i"p_L) elR" relative to ^1, ^, when /ipy 

is identified with hpc and /ip _l is identified with hp_[_ = hpc±- In a similar 
form, we decompose vertical SO{m +1)/ SO{m) matrices 



vp3 



— vp 



wp 



+ 



(wpiht) 
-(^vp||,0")^ 

(0,^;^^) 
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into tangential and normal parts relative to e^x via corresponding decom- 
positions of h-vectors vp = {vp^.v'pj^) gM™ relative to ^1, ^, when vp\\ is 

identified with vpc and t; p j_ is identified with vp± = vpc-^- 

There are analogous matrix decompositions relative to the flow direction. 
In the h-direction, we parametrize 



(/ley, h~e±) 



— {he\\, h'e±) hO 
when e/iY e hp, [he\\,h'e±) e W and h~e± e and 



(0, 
-(0,'0)^ huir 



e so{n + 1), 



where 



hW-r — 



w 

-w^ he 



e so{n), w e he eso{n - 1). 



In the v-direction, we parametrize 



(ve\\,v*e±) 



when e^Y G "^P, {ve\\,v^ ±) G M™ and v 



e I G 



pm—l 



, and 



where 



r^Y=7i)Yjc 



w 



0^ (0, ) 

-(0, O")^ VWr 



G so(m + 1), 



G S0(m), G R"" , f0 G50(m - 1). 



(31) 



The components he\\ and h~e± correspond to the decomposition 

ehY = /ig(7^,7/)efcx + (7r)^J^e± 

into tangential and normal parts relative to e^x- In a similar form, one 
considers ven and v e ± corresponding to the decomposition 



e^Y = 'yg(7^,7;)e^x + {lr)^\^^^- 
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Using the above stated matrix parametrizations, we get 



he± 

[r/iY, G/iy] 

and 

ve± 




he± 





h~e 



with 



(32) 







{0,w 




(^—~v-h'e±,he\i'v') 
(^—lj'-h~e±,he\\l^) hO 



G hp; 





fe_L 



iv~ej_ 




Gso(m + l), with 
G so(m); 



(33) 



f e _|_ 







G vp±; 



(— V ■ f V_L, i;e|| V) 

— (— V ■ f V_L, t>e|| V) 



G t;p. 



We can use formulas flHzl) and flHHl) in order to write the structure equa- 
tions ( 124|) and ( l25i) in terms of N-adapted frames and connection 1-forms 
soldered to the Klein geometry of N-anholonomic spaces using the relations 
(130|1 . One obtains respectively the G-invariant N-adapted torsion and cur- 
vature generated by the canonical d-connection, 

T(7,, 70 = (Dx7r - Dy7/)J e = DxBy - DyBx + [Fx, by] - [Fy, ex] (34) 
and 



R(7,, 70e = [Dx, Dy] e = DxFy - DyFx + [Fx, 



(35) 



where ex == 7je, Gy == 7rje, Fx == 7ijr and Fy == 7Tjr. The formu- 
las ( IMI) and ( l35l) are respectively equivalent to ( |T3|) and ( 1721) . In general, 
T(7r, 7z) 7^ and R(7r, 7;)^ cannot be defined to have constant matrix coef- 
ficients with respect to a N-adapted basis. For N-anholonomic spaces with 
dimensions n = m, we have T(7t-, 7^) = and R(7t, 7z)e defined by constant, 
or vanishing, d-curvature coefficients (see discussions related to formulas (17^ 
and ( |70l) ). For such cases, we can consider the h- and v-components of ( l34l) 
and (1351) in a similar manner as for symmetric Riemannian spaces but with 
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the canonical d-connection instead of the Levi Civita one. One obtains, 
respectively, 



(D;,x7r-D^Y70J^e (36) 



and 



/iR(7^, 7i)/;,e = [Bhx, ^hy] he = Bhx^'hY - D/iyL^x + [L/iX, L^y] (37) 
wR(7^, -fi)ve = [D^x, D^y] ^^e = Di,xCi,y - D^yC^jx + [C^x, C^y] ■ 

Following N-adapted curve flow parametrizations ( 132|) and ( l33i) , the equa- 
tions and (!H7|) are written 

= D/ix^e|| + 1? ■ /i"e'_L, = D^x'^en + V ■ f ; (38) 
= w — henlj' + 'Dhxh'e ±, = w — ven^ + D^xv^e ±; 



and 
D 



'hx^ - D/tY"i^ + Z?-© = h~e ±, Dyx.^ - 'DyY^v + *v \v& =v*e ±; 
^hxh@—~u + = 0, Dj,xf ©— V (8)& + ro(8)V = 0. (39) 

The tensor products and interior products are defined in the form: for the 
h-components, ® denotes the outer product of pairs of vectors (1 x n row 
matrices), producing n x n matrices A <^ B = B , and J denotes mul- 
tiplication of X matrices on vectors (1 x n row matrices), such that 

A\ ^i?®cj = (^A ■ B^ C which is the transpose of the standard matrix 
product on column vectors. Likewise, for the v-components, we just change 
n — i> m and A ^ A. For the sequel, A will denote the skew product of 
vectors AAB = A(g)B-B(g)A. 

The variables ey and 0, written in h- and v-components, can be ex- 
pressed in terms of the variables 1^ ,w , h e ± and V, ro, f e ^ (see respec- 
tively the first two equations in (1381) and the last two equations in ( l39ll ) 

by 

/ley = -D^x(^ ■ h~e±)^ ^e^ = -D^-^C^ " ^ e"^), 

and 

h& = ® — (g) ~u) , vG = D,^x ® ^ — ^ ® V) . 
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Substituting these expressions, correspondingly, in the last two equations in 
and in the first two equations in flHIH) . we express 



contained in the h- and v-flow equations respectively on Ij' and V, consid- 
ered as evolution equations via DhY^ = and D^y^ = 

~Ur = 'Dhx.'^ — ^\'Df^^{l^ ®w — w ®~u) — Rh~ej_, (40) 
*Vr = 'Dvx'^ — '*^jD~x ®w — ^ ® V) — S v*e_i, 

where the scalar curvatures of the canonical d-connection, R and S are 
defined by formulas (17^ in Appendix. For symmetric Riemannian spaces 
like SO{n + 1)/ SO{n) ~ S*"", R is just the scalar curvature x = 1; see [T7] . 
On N-anholonomic manifolds under consideration, the values R and S are 
constrained by be certain zero or nonzero constants. 

The above presented considerations lead to a proof of the following main 
results. 

Lemma 4.1 On N-anholonomic spaces with constant curvature matrix co- 
efficients for the canonical d-connection, there are N-adapted Hamiltonian 
symplectic operators, 

hj = T>h^ + B^^{ir-)^ and t;:^ = D,x + D;^ (V-) V, (41) 

and cosymplectic operators 

= D;,x + ^jD^i(^A) and = D^x + VjD;^ (Va) . (42) 

The properties of operators fBTj) and fH21) are defined by 

Theorem 4.1 The d-operators J = {hj',vj') and Ti = {HTi^vl-C) are res- 
pectively {0{n — 1), 0{m — l))-invariant Hamiltonian symplectic and cosym- 
plectic operators with respect to the flow variables (^, V). Consequently the 
curve flow equations on N-anholonomic manifolds with constant d-connection 
curvature have a Hamiltonian form: the h-flows are given by 

'v'r = hTi {w) — R h~e j_ = h9\ [h~e ±) — R h~e±, 
w = hj (/i"e x) ; (43) 

the v-fiows are given by 

*Vr = vH (w) — S v*e ± = vDl (v^±) — S v^±, 
w = vj e"_L) , (44) 



27 



where the hereditary recursion d-operator has the respective h- and v-compo- 
nents 

h^ = HHo hj and = vH o vj. (45) 

Proof. Details in the case of holonomic structures are given in Ref. [13] 
and generalized to arbitrary gauge groups in ^34j- For the present case, the 
main additional considerations consist of soldering certain classes of general- 
ized Lagrange spaces with {0{n — l),0(m — l))-gauge symmetry onto the 
Klein geometry of N-anholonomic spaces. □ 

4.2 bi-Hamiltonian anholonomic curve flows and soli- 
tonic hierarchies 

Following the usual solitonic techniques, see details in Refs. [171 [35], the 
recursion h-operator from fl45p . 

km = D;,x(D/,x + D^i(l?-)^) +iy^jD^^(^AD^x) (46) 
= + |D,x|' + (IT-) - ^jD,i(^,A), 

generates a horizontal hierarchy of commuting Hamiltonian vector fields 
h~e^^^ starting from h~e^f^ = 'v'l given by the infinitesimal generator of 
/-translations in terms of arclength I along the curve. 

A vertical hierarchy of commuting vector fields v*e^^ starting from v*e'^^ 
= is generated by the recursion v-operator 

vm = d„x(d,x + d;^(V-)V) + Vjd;^(Vad,x) (47) 

= + |D,x|' + (V-) V; - VjD;i(V,A). 

There are related hierarchies, generated by adjoint operators 9^* = {h^i*, 
v^*), of involutive variational h-covector fields w^'^^ = 5 {hH^'^'>^ /5'v' in 
terms of Hamiltonians hH = hH^''\~i^,lTi,lj'2h-- ) starting from ro'-^^ = 
~v,hH^^^ = and of involutive variational v-covector fields w'^''^ = 

5 (fiJ'^'^)) / 5V in terms of Hamiltonians vH = vH^''\^, V;, V2/, . . .) start- 
ing from w^^^ = *v',vH^^^ = The relations between hierarchies are 
established correspondingly by formulas 

/i-i^f = hn (^('^^ ^('=+1)) = hj (h-^^l^^ 

and 
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where k = 0, 1,2, . . .. All hierarchies (horizontal, vertical and their adjoint 
ones) have a typical mKdV scaling symmetry, for instance, I—* XI and ^ 
A~^^ under which the values h~e^^^ and hH^^^ have scaling weight 2 + 2 A;, 
while w^'^'^ has scaling weight 1 + 2k. 

The above presented considerations prove 

Corollary 4.1 The recursion d-operator ( fT^j j gives rise to N-adapted hi- 
erarchies of distinguished horizontal and vertical commuting bi-Hamiltonian 
flows, respectively, on ^ and V that are given by 0{n — 1)- and 0{m — 1) 
-invariant d-vector evolution equations, 

with horizontal Hamiltonians /li/^'^+i'^) = /i/^C^'+i'^) — R hH^^'^^ and 

with vertical Hamiltonians vH^^^'^'^'^ = vH'^^^^'^^ — S vH^^'^\ for k = 
0,1,2,.... The d-operators Ti and J are N-adapted and mutually com- 
patible, from which one can construct an alternative (explicit) Hamiltonian 
d-operator Q = H o J' oTi =£Ho7i. 

4.2.1 Formulation of the main theorem 

The main goal of this paper is to prove that for any regular Lagrange sys- 
tem we can define naturally a N-adapted bi-Hamiltonian hierarchy of flows 
inducing anholonomic solitonic configurations. 

Theorem 4.2 For any anholonomic vector bundle with prescribed d-metric 
structure, there is a hierarchy of bi-Hamiltonian N-adapted flows of curves 
'j{t,1) = [h'~f{T,l),v'~f{T,l)] described by geometric nonholonomic map equa- 
tions. flows are defined as convective (traveling wave) maps 

(h-f)^ = {h-f)^y^ and (177)^ = (v-f)^-^ . (48) 

The flows are defined by non-stretching mKdV maps 

-M. = DL(M.x + ^|D/.x(/^7),xIL (Mhx. (49) 
-{v-i)^ = D2x(t;7)^x + 2 l°^x(t^7)^,xl^g M^,x > 
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and the +2, . . . flows as higher order analogs. Finally, there are —1 flows 
are defined by the kernels of recursion operators ( [7^ and inducing non- 
stretching wave maps 

DhY (/i7)/.x = and D„y (t^7).x = 0- (50) 
We carry out the proof in the next section I4.2.2I 

Remark 4.1 Counterparts of N- adapted bi-Hamiltonian hierarchies and re- 
lated solitonic equations in Theorem \4.2\ can be constructed for SU{n) © 
SU{m) / SO{n) © SO{m) as done in Ref. JJ^ for Riemannian symmetric 
spaces SU{n)/ SO{n). Such results may be very relevant in modern quantum 
/ (non) commutative gravity. 

Indeed, constructions similar to Theorem 14.21 can be carried out in gravity 
models with nontrivial torsion and nonholonomic structure and related ge- 
ometry of noncommutative/ superspaces and anholonomic spinors. Finally, 
it should be emphasized that a number of exact solutions in gravity [251 ES] 
can be nonholonomically deformed in order to generate nonholonomic hi- 
erarchies of gravitational solitons of type ( 148|) . ( 149|) or ( l50i) . which will be 
considered in further publications. 

4.2.2 Proof of the main theorem 

We provide a proof of Theorem 14.21 for the horizontal flows. The approach 
is based on the method provided in Ref. [TH |T7j, but in this proof the Levi 
Civita connection on symmetric Riemannian spaces is replaced by the hori- 
zontal components of the canonical d-connection in a generalized Lagrange 
space with constant d-curvature coefficients. The vertical constructions are 
similar. 

One obtains a vector mKdV equation up to a convective term (which can 
be absorbed by redefinition of coordinates) defining the +1 fiow for /i e _l = 

when the +(fc + 1) fiow gives a vector mKdV equation of higher order 3 + 2fc 
on ^ and there is a h-fiow = l^i arising from /i e ^ = and he \\ = 1 
belonging outside the hierarchy generated by h^K. Such fiows correspond to 
N-adapted horizontal motions of the curve 7(r, /) = [/i7(r, Z), f7(r, I)] given 

by 

(M. = / {{hl)^^ , D,x (/^7),x ' DL (M,x - ■ •) 
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subject to the non-stretching condition | {h'y)^^-^ \fig = 1, where the equation 
of motion is derived from the identifications 

and so on, which maps the constructions from the tangent space of the curve 
to the space hp. For such identifications we have 



[L/ix, [Lfex, G/ix]] 





(0,^ 



(0,^ 
hO 



v\ , 



G hp, 
hO 



and so on, see similar calculus in ( !32l) . At the next step, noting for the +1 
h-flow that 



h~e j_ = ~ui and h~e\\ = — D^^ ' '^i) = ^2 



we compute 



{he\\, h~e 



(/ie||, h~ej_^ 



hO 



- 0, 



1, 



hO 



+ 







- (0, irhx.f hO 



hXJ 



— D/ix [L/tx, G/ix] + - [L/ix, [L/iX7 G/ix]] 

3 

= -^^h.x [Lfex, dhx] - -l^l^e/ix- 

From the above presented identifications related to the first and second terms, 
it follows that 

I w"!^ = < [L/iX, e/,x] , [L/ix, e^x] >hp< — > /^g (D/,x (/i7)hx > D^x (^7)hx) 
= |Dhx(/i7)hxlhg' 
and we can identify Vhx. [L/iX, e/ix] with D^^ (^7)/iX ^^"^ write 



3 

-GhY < — > (/i7);,x + 2 (Mhxihg ih^)hx 
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which is just the first equation in the Theorem 14.21 defining a non- 
stretching mKdV map h-equation induced by the h-part of the canonical 
d-connection. 

Using the adjoint representation ad (■) acting in the Lie algebra hQ = 
hp © so(n), with 



arf([L/jx,e/ix])e/ix 





o,"o 



T 



0, 



G so{n + 1), 



where 







-v'^ hO 

we note (applying ad{[Lhx,ehx]) again ) 



G 3o{n) 



a(i([L/,x,e/ix])^e/iX = -I^T 







- 1, 



1, 







- V e/,x- 



Hence equation fl49p can be represented in alternative form 



- (^7). = (/^7),x - ^R-'ad (D,x {hj),^? (/^7),x , 

which is more convenient for analysis of higher order fiows on Ij' subjected 
to higher-order geometric partial differential equations. Here we note that 
the flow on Ij' corresponds to just a convective (linear traveling h-wave 
but subjected to certain nonholonomic constraints ) map equation fHHj) . 

Now we consider a —1 flow contained in the h-hierarchy derived from 
the property that h e ± is annihilated by the h-operator hj^ and mapped 
into hy\{h~e ±) = 0. This mean that hj'{h~e±) = w = 0. Such properties 
together with (13T|) and equations (140|) imply = and hence hV^ehy^ = 
[L^, e/ix] = for hV^ = hD^ + [Lt, ■]• We obtain the equation of motion for 

the h-component of curve, h'y{T,l), following the correspondences DhY * ^ 

hVr and h'ji < — > e/^x, 

D,,Y(/i7(r,/)) = 0, 
which is just the flrst equation in ( l50l) . 



Finally, we note that the formulas for the v-components, stated by The- 
orem 14.21 can be derived in a similar form by respective substitution in the 
above proof of the h-operators and h-variables into v-ones, for instance, 
h'j vy, h~e± — > v*e±, ^ V, ^ — > w,D/jx — > D,;x, D^y 
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4.3 Nonholonomic mKdV and SG hierarchies 



We present explicit constructions when solitonic hierarchies are derived fol- 
lowing the conditions of Theorem I4.2[ 

The h-flow and v-flow equations resulting from ( l50i) are 

l^r = —Rh~e± and Vr = — 5'fVj_, (51) 

when, respectively, 

= = —T^hxh~e j_ + he^^~v, D/jx/iey = h~e j_ ■ 'v' 

and 

= = — D^jx'w^x + fe|| V, D^jx'weii = v% x_ ■ V. 
The d-flow equations possess horizontal and vertical conservation laws 

Dftx((/ie,|)2 + |/i"?xr) =0, 
for (/ie||)2 + \h~e±\'^ =< her, >hp= \ {h'y)^ and 

D,Y((t^e||)2 + |i;V^|2) =0, 
for {ve^iY + =< veT-,veT >^p= | (v'j)^ This corresponds to 

Dhxl (/i7)r = D^xl (f7)r l^g = 0- 

(The problem of formulating conservation laws on N-anholonomic spaces — 
in particular, on nonholonomic vector bundles — is analyzed in Ref . [26] . In 
general, such laws are more sophisticated than those on (semi) Riemannian 
spaces because of nonholonomic constraints resulting in non-symmetric Ricci 
tensors and different types of identities. But for the geometries modeled for 
dimensions n = m with canonical d-connections, we get similar h- and v- 
components of the conservation law equations as on symmetric Riemannian 
spaces.) 

Without loss of generality we can rescale conformally the variable r in 
order to get | (hj)^ \\^ = 1 and (perhaps by a different rescaling) | (^7)^ = 
1, so that 

(/ie||)^ + \h~e±\'^ = 1 and (fey)^ + jt'^xl^ = 1. 

In this case, we can express /ley and h~e± in terms of 1? and its derivatives 
and, similarly, we can express fey and v^± in terms of V and its derivatives, 
which follows from (15T1) . The N-adapted wave map equations describing the 
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— 1 flows reduce to a system of two independent nonlocal evolution equations 
for the h- and v-components, 



r = -Dri K/i?2 _ 1- 12 - and V. = -D"i V S 



For N-anholonomic spaces of constant scalar d-curvatures, we can rescale the 
equations on r to the case when the terms R^, S"^ = 1, and the evolution 
equations transform into a system of hyperbolic d-vector equations, 



'Dh:s.{v r) = -y^ - \v tI"^ V and D^x( r) = -y 1 - I , (52) 

where Dhx = dhi and D^x = d^i are usual partial derivatives, respectively, 
along hi and vl with v r and v r considered as scalar functions for the 
covariant derivatives D/^x and Di,x deflned by the canonical d-connection. 
It also follows that he and v e _[_ obey corresponding vector sine-Gordon 
(SG) equations 



{l-\h-§^\')~^ dhi{h-§i_)] =-h-^± (53) 



and 

(^y/(l-|i;V^|2)-i5,i(t;V^)^ = -yV^. (54) 

The above presented formulas and Corollary 14. II imply 

Theorem 4.3 The recursion d-operator 9^ = (/i9^,/i9^) generates two 
hierarchies of vector mKdV symmetries: the first one is horizontal (see ^^), 

l/f = 'vhu^^r'^ = h^{-Vhi) = l^m + \\'v\^l^hu (55) 



V 



(2) 



5 



5 / / 1 — y 1 9\ I — >■ I 2 ^1 — > 1 4 \ — > 1 I — > I 2 — ^ 

+ 2 ( (I ^ I + I ^^il +4!^ I ) Vhl--\Vhl\ V, 

and so on, 

with all such fl,ows commuting with the —1 fl^ow 

{^r)~^ = h~e± (56) 
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associated to the vector SG equation [53^) : the second one is vertical (see 

3 

5 

and so on, 

with all such flows commuting with the —1 flow 

= v"^^ (58) 

associated to the vector SG equation ( fJ^p . 

Furthermore, the adjoint d-operator 9^* = J oTi generates a horizontal 
hierarchy of Hamiltonians, 

2' ' 2' ^" 8 



/,^(2) _ ^Tt 12 -,T-|2 



2I ^ 2M| -^1^1 \Vhl\ - ■ V hl} + Z^\V 



2 1 2m| - 
and so on, 

and vertical hierarchy of Hamiltonians 

2' ' ' 2' ' 8 

1 I 4^ I o 3 I 4 I o I 4 I o 1 / 4 4 ^ 1 

2l^2^«| - 
and so on 



vH^'^ = t;i7« = --|V.n' + ^|Vr, (60) 

Z Z o 



all of which are conserved densities for respective horizontal and vertical 
— 1 flows and determining higher conservation laws for the corresponding 
hyperbolic equations fl^ and fl3^ . 

The above presented horizontal equations fl3H|) . fl33|) . fl3U|) and fl3IH) and 
of vertical equations ( 15^ . ( !57|) . (l58l) and ( !60|) have similar mKdV scaling 
symmetries but with distinct parameters Xh and because, in general, there 
are two independent scalar curvatures R and S , see ( I7il) . The horizontal 
scaling symmetries are hl^Xhhl~v — > {\h)~^ ~v and r — > (A/i)^'*'^'^ , for A; = 
— 1, 0, 1, 2, and likewise for the vertical scaling symmetries. 
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Finally, we consider again the Remark 14.11 stating that similar results 
(proved in Section 4) can be derived for unitary groups with complex vari- 
ables. The generated bi-Hamiltonian horizontal and vertical hierarchies and 
solitonic equations will be different from those defined for real orthogonal 
groups; for holonomic spaces this is demonstrated in Ref. [H]. This distin- 
guishes substantially the models of gauge gravity with structure groups like 
the unitary one from those with orthogonal groups. 

5 Applications in Geometric Mechanics, Fin- 
sler Geometry and Gravity 

Here we consider some interesting examples when the data defining funda- 
mental geometric structures in mechanics and Finsler geometry, or exact so- 
lutions in gravity, can be transformed into solitonic hierarchies. The possibil- 
ity to describe gravitational and electromagnetic interactions by mechanical 
models, and vice versa, will be investigated. 

5.1 Geometric mechanics and Finsler geometry 

A differentiable Lagrangian L{x,y), i.e. a fundamental Lagrange function, 
is defined by a map L : {x,y) G TM — >■ L{x,y) G M of class C°° on TM = 
TM\{0} and continuous on the null section : M — > TM of vr. A regular 
Lagrangian has non-degenerate Hessian 

when rank \gij\ = n on TM. 

Definition 5.1 A Lagrange space is a pair L" = [M, L{x,y)] with ^^^gij 
being of fixed signature over V = TM. 

The notion of Lagrange space was introduced by J. Kern [1] and elab- 
orated in details by R. Miron's school, see Refs. [21 El EZ], as a natural 
extension of Finsler geometry [3HI ESI EHl 123] (see also Refs. [2ni SO], on 
Lagrange-Finsler super/noncommutative geometry). Straightforward calcu- 
lations (on nonholonomic manifolds they are reviewed in Refs. [26]) establish 
the following results: 
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1. The Euler-Lagrange equations 

d /dL\ dL _^ 

dr \dy^ J dx^ 

where ?/* = ^ for x*(r) depending on parameter r, are equivalent to 
the "nonhnear" geodesic equations 

defining paths of a canonical semispray 

■ d d 

where 

2G'(x, = - (^V^Y-^y'^ - 1^ 
^ ^' 2 ^ \dy'dx^^ dx' 

with ^^^g"^^ being inverse to flUTl) . 

2. There exists on V ~ TM a canonical N-connection 



^L)^^ ^ dG\x,y) 

3 Qyi 



(62) 



defined by the fundamental Lagrange function L(x, y), which prescribes 
nonholonomic frame structures of type (jlj) and ([5]), '^^^e^ = (ej, Ca) and 
(L)e/. = (e\e''). 

3. There is a canonical metric structure 

(^)g = Qijix, y) e' ® + g^j{x, y) ® e^' (63) 
constructed as a Sasaki type lift from M for gij{x, y). 

4. There is a unique metrical and, in this case, torsionless canonical d- 
connection (^^D = {hD,vD) with the nontrivial coefficients with re- 
spect to ^^^e^ and (^)e^ parametrized respectively F"^^ = {D-j^,C^^, 
for 

L'jk = ^9'''{Gkgjh + ejgkh-ehgjk), (64) 
1 



C'jk = i(^k9jh + ebgkh - eegbc) 

defining the generalized Christoffel symbols, where (for simplicity, we 
omitted the left up labels (L) for N-adapted bases). The connection 
•^^^D is metric compatible and torsionless, see an outline of main defi- 
nitions and formulas in the next subsection and Appendix. 
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We conclude that any regular Lagrange mechanics can be geometrized as 
a nonholonomic Riemann manifold V equipped with canonical N-connection 
( !62l) and adapted d-connection (l6^ and d-metric structures ( l63|) all in- 
duced by a L{x,y). In some approaches to Finsler geometry and general- 
izations [32], one consider nontrivial non- metric structures. For instance, 
the so called Chern d-connection '-'^D is also a minimal extension of the 
Levi Civita connection when h '-^'^D((yf) =0 but v '-^'^D(/i) 7^0, see formulas 
(fT6l) . Such generalized Riemann-Finsler spaces can be modeled on nonholo- 
nomic metric-affine manifolds, with nonmetricity, see detailed discussions in 
[261 [2], [3], [21] . We note that the N-connection is induced by the semispray 
configurations subjected to generalized nonlinear geodesic equations equiva- 
lent to the Euler-Lagrange equations. A N-connection structure transforms 
a Riemannian space into a nonholonomic one with preferred honholonomic 
frame structure of type (jlj) and ([5]). 

Remark 5.1 Any Finsler geometry with a fundamental Finsler function 
F{x,y), being homogeneous of type F{x, Xy) = XF{x,y), for nonzero A G M, 
may be considered as a particular case of Lagrange geometry when L = F^. 
We shall apply the methods of Finsler geometry in this work^ 

For applications in optics of nonhomogeneous media and gravity (see, 
for instance, Refs. [31 [2Sl IHl ES]) considers metric forms of type gij ~ 
gA(a:,j/) Wg-j(^x,y) which cau not be derived explicitly from a mechanical 
Lagrangian. In the so-called generalized Lagrange geometry one considers 
Sasaki type metrics ( l63l) with certain general coefficients both for the d- 
metric and N-connection and canonical d-connection, i.e. when ^^^gij — >■ 
[gij(x,y), hab{x,y)], and ^^^iVj — > Nj{x,y). We shall use the term (gener- 
alized) Lagrange-Finsler geometry for all such geometries modeled on tan- 
gent bundles or on arbitrary N-anholonomic manifold V. In original form, 
such spaces were called generalized Lagrange spaces and denoted GL" = 
{M,gij{x,y)), see [21 E]. 

5.2 Analogous Models and Lagrange-Finsler geometry 

For a regular Lagrangian L{x,y) and corresponding Euler-Lagrange equa- 
tions it is possible to construct canonically a generalized Finsler geometry on 

^In another direction, there is a proof |32] that any Lagrange fundamental function 
L can be modeled as a singular case in a certain class of Finsler geometries of extra 
dimension. Nevertheless the concept of Lagrangian is a very important geometrical and 
physical one and we shall distinguish the cases when we model a Lagrange or a Finsler 
geometry: A physical or mechanical model with a Lagrangian is not only a "singular" case 
for a Finsler geometry but reflects a proper set of geometric objects and structures with 
possible new concepts in physical theories. 
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tangent bundle [Il[21[3j. In our approach, such constructions follow from the 
Theorem 12. II if instead of Qajsix) there are considered the metric components 



1 d^L{x, y) 

2 dy-dyi- ' 



(65) 



introduced in ([T]) for e^{x,y) = 5^. For arbitrary given gij{x,y) and 
e^-(x, y) on TM, defining a d-metric ffT7|) . we can compute the corresponding 
N-connection ([H]) and canonical d-connection coefficients flUUj) . In general, 
the coefficients of d-curvatures fl7^ are not constant and we are not able 
to generate solitonic hierarchies. In order to solve the problem we have to 
choose e^-(x, y) so that the equations 



are satisfied (similarly to (|T8l) . for given L and Qab), resulting in constant 
curvature coefficients with respect to N-adapted bases. This allows us to 
define D)^x fro^i formulas (jH]), fH^ or fl^Ul) . as pairs of h- and 

v-operators with "nonholonomic mixture" (cf. Theorem 14.21) . 

Corollary 5.1 Any regular Lagrangian \-{x,y) induces a hierarchy of bi- 
Hamiltonian N-adapted flows of curves described by geometric nonholonomic 
(solitonic) map equations. 

We can describe equivalently any regular Lagrange mechanics in terms of 
solitonic equations, their solutions and symmetries. 

Remark 5.2 Considering L{x,y) = F'^{x,y), where the homogeneous 
F{x,Xy) = XF{x,y) is a fundamental Finsler metric function (see, for in- 
stance, J^), we can encode in solitonic constructions all data for a Finsler 
geometry. 

Here it should be noted that there were elaborated alternative approaches 
when gravitational effects are modeled in continuous and discrete media, see 
review [M], but they do not allow one to generate the Einstein equations 
starting from a Lagrangian in a mechanical models like in [H [21 E], or from a 
distribution on a nonholonomic manifold [301 [S], EH |26] . I^i our approach, for 
a d-metric ( |T71) induced by ( l65i) (introduced in ([1])), we can prove a geometric 
mechanical analog of Theorem 12.31 

Theorem 5.1 Any regular Lagrangian L{x,y) and frame structure e^{x,y) 
define a nonholonomic (semi) Riemannian geometry on TM. 
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This way we can model various gravitational effects by certain mechanical 
configurations. 



Corollary 5.2 A regular Lagrangian L(x, ?/) and frame structure e^{x,y) 
on TM model a vacuum gravity configuration with effective metric cjap = 
[9ijy9ij] p7| j if the corresponding canonical d-connection ^UM) has vanishing 
Ricci d-tensor [7^ . 

We conclude that for a fixed Lagrangian L{x,y) we can define frame 
structures e„-(x, y) on TM inducing Einstein spaces or curved spaces with 
constant curvature coefficients (computed with respect to certain N-adapted 
basis). 

Definition 5.2 A geometric model defined by data E"^"^ = [gij{x'^ ,y'^) , 
Qabix^i y^)-, N°'{x^, y'^)] is N-anholonomically equivalent to another one given 
by data E""*"™ = lgij{x'' ,y'^),gab{x'' ,y'^), N^{x^ ,y'^)] if for the same splitting 
n + m withea — > Cq, there are nontrivial polarizations [riij{x^ ,y'^),riab{x'^ ,y'^) , 
7]'^{x^,y'')] for which g^p = Vafsdap and = r^fiVf . 

In general, the physical and geometric properties of two such N-anholono- 
mically related spaces E""*""^ and E""*""^ are very different. Nevertheless, we 
are able to compute the coefficients of geometrical and physical objects, define 
symmetries, conservation laws and fundamental equations on a space from 
similar values of another one if the polarizations are stated in explicit form. 
In some particular cases, we can consider that E""*""^ is defined by the data 
with constant curvature coefficients with respect to an N-adapted basis (or 
an exact solution of the Einstein equations) but E"-"*""^ is related to a regular 
Lagrange/ Finsler geometry. 

5.3 Modeling field interactions in Lagrange— Finsler ge- 
ometry 

Let us consider a regular Lagrangian 



where mo, cq are constants, A^{x) is a vector field and a^j{x) = a.ij{x) is a sec- 
ond rank symmetric tensor. The metric coefficients f l65l) are gij = mo dijix) 
and their frame transforms ([1]) are given by 



l{x,y) = mQaJx)y'y^ + eQAi{x)y\ 



(66) 



9ab{x,y) = e^-{x,y) ef{x,y) gab{x) = niQ aab{x,y), 



(67) 
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where aab = e^- e^- agb. For simplicity, we consider frame transforms to 
local bases Cq = (ei,ea) when = e^- d/dy- are holonomic vectors but 
Cj = Cj- d/dx- will be defined by a N-connection like in (jl]). We compute 
(see Theorem 12.11 and formula ([H])) 

& = l 'r,,y'y' + —a^' ^F,, y\ 

where jk are Christoffel symbols of the tensor gab constructed by using 
derivatives d/dx^, ^Fj^ = ^( ^Fi^Aj — ^DjA^) with the covariant derivative 
3Dj defined by ^r.j^, and 

N'^ix,y)= 'r^,y'- 3F,k. (68) 
The N-connection curvature ^ is 

"jfc — y ^ bjk ^[k ^ i], 

where ^R^j^ is the curvature tensor of the Levi Civita connection for gqb{x), 
computed with respect to a usual coordinate base and then nonholonomically 
transformed by e^{x,y). 

For mo = mc and Cq = 2e/m, where m and e are respectively the mass 
and electric charge of a point particle and c is the light speed, the Lagrangian 
(1661) describes the dynamics of a point electrically charged particle in a curved 
background space M with metric ajj(x). The formula ( |671) states a class of 
nonholonomic deformations of the metric on M to TM. We can treat ^Fj^. 
as an effective electro-magnetic tensor field on TM. The auto-parallel curves 
u"{t) = [x*(r), |/*(r) = (ix'^/r], parametrized by a scalar variable r, adapted 
to N'j{x,y) (EHD, are 

^+ ^T\^{x,y)yV= 'F\{x,y)y\ 

The nonholonomic Riemannian mechanical model of L(x, ?/) is completely 
defined by a d-metric ( |T7I) with coefficients f l67I) . 

g = mo ajj(x, e* ® e-^ + mo aij{x, y) e* ® e\ 

where e* are elongated by N^j from ( l68l) . From Corollary 15.21 we have: 

Conclusion 5.1 Any solution of the Einstein-Maxwell equations given by 
gravitational field a^j{x) and electromagnetic potential A^{x) is N-anholono- 
mically equivalent to a mechanical model with regular Lagrangian and 
associated nonholonomic geometry on TM with induced d-metric g and N- 
connection structure N^j. 
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If the frame coefficients e^- are chosen to yield constant d-curvatures 
coefficients, we get a particular statement of the main Theorem I4.2t 

Corollary 5.3 The Einstein-Maxwell gravity equations contain solitonic so- 
lutions that can be modeled as a hierarchy of bi-Hamiltonian N-adapted flows 
of curves. 

In a similar form, we can derive spinor-solitonic nonholonomic hierarchies 
for Einstein-Dirac equations if we consider unitary groups and distinguished 
ChfFord structures [2S]- 

6 Conclusion 

In this paper, it has been shown that the geometry of regular Lagrange me- 
chanics and generalized Lagrange-Finsler spaces can be encoded in nonholo- 
nomic hierarchies of bi-Hamiltonian structures and related solitonic equations 
derived for curve flows. Although the constructions are performed in explicit 
form for a special class of generalized Lagrange spaces with constant ma- 
trix curvature coefficients computed with respect to canonical nonholonomic 
frames (induced by Lagrange or Finsler metric fundamental functions), the 
importance of the resulting geometric nonlinear analysis of such physical sys- 
tems is beyond doubt for further applications in classical and quantum field 
theory. This remarkable generality appears naturally for all types of models 
of gravitational, gauge and spinor interactions, geometrized in terms of vec- 
tor/ spinor bundles, in supersymmetric and/or (non) commutative variants 
when nonholonomic frames and generalized linear and nonlinear connections 
are introduced into consideration. 

If the gauge group structure of vector bundles and nonholonomic mani- 
folds is defined by an orthogonal group H = SO{n) © SO{m) acting on the 
base/ horizontal and typical fiber /vertical subspaces with, say, respective 
dimensions of n and m, then these subspaces will be Riemannian symmetric 
spaces G/H whose geometric properties are determined by the Lie groups 
G = SO{n) © SO{m) and G = SU{n) © SU{m). This structure can be for- 
mulated equivalently in terms of geometric objects defined on pairs of Klein 
geometries. The bi-Hamiltonian solitonic hierarchies are generated naturally 
by recursion operators associated with the horizontal and vertical fiows of 
curves on such spaces, with these fiows (of SG and mKdV type) being geo- 
metrically described by wave maps and mKdV analogs of Schrodinger maps. 

In the case of holonomic manifolds, the construction of bi-Hamiltonian 
curve fiows, soliton equations, and geometric PDE maps has recently been 
generalized in Ref. [34i to all Riemannian symmetric manifolds M = G/H 
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including all compact semisimple Lie groups K = G/H as given by G = 
K X K and H =diagG'. This generalization makes it possible to study more 
general field models, such as quaternionic and octonian systems. 

Considering arbitrary curved spaces and mechanical or field systems, it 
is not clear if any bi-Hamiltonian structures and solitonic equations can be 
derived from curve fiows. The answer seems to be negative because arbitrary 
curved spaces do not possess constant matrix curvatures with respect to cer- 
tain orthonormalized frames which is crucial for encoding recursion operators 
and associated solitonic hierarchies^ Nevertheless, the results proved in our 
work provide a new geometric method of solitonic encoding of data for quite 
general types of curved spaces and nonlinear physical theories. This follows 
from the fact that we can always nonholonomically deform a curved space 
having nonconstant curvature given by the distinguished linear connection 
into a similar space for which this curvature becomes constant. 

Our approach employs certain methods elaborated in Finsler and La- 
grange geometry when the geometric objects are adapted to nonholonomic 
distributions on vector/ tangent bundles, or (in general) on nonholonomic 
manifolds. This accounts for existence of the canonical nonlinear connection 
(defined as non-integral distributions into conventional horizontal and verti- 
cal subspaces, and associated nonholonomic frames), metric and linear con- 
nection structures all derived from a regular Lagrange (in particular, Finsler) 
metric function and/ or by moving frames and generically off-diagonal met- 
rics in gravity theories. Of course, such spaces are not generally defined to 
have constant linear connection curvatures and vanishing torsions to which 
the curve-fiow solitonic generation techniques can not applied. 

We proved that for a very large class of regular Lagrangians and so- 
called generalized Lagrange space metrics, there can be nonholonomic de- 
formations of geometric objects to equivalent ones on generalized Lagrange 
configurations, for instance, with constant Hessian for the so-called absolute 
energy function. The curvature matrix, with respect to the correspondingly 
adapted (to the canonical connection structure) frames, can be defined with 
constant, even vanishing, coefficients. For such configurations, we can apply 
the former methods elaborated for symmetric Riemannian spaces in order to 
generate curve-fiow solitonic hierarchies. Here should be noted that such spe- 
cial classes of nonholonomic manifolds are equipped with generic off-diagonal 
metrics and possess nontrivial curvature for the Levi Civita connection. They 
are characterized additionally by nontrivial nonlinear connection curvature 

^"a series of recent works [HI |331 Hll SB ESI SZ] provide a number of generalizations and 
examples when solitonic hierarchies can be derived for arbitrary Einstein, Riemann-Cartan 
and Larange-Finsler spaces, their Ricci flow evolutions and/or quantum deformations 
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and nonzero anholonomy coefficients for preferred frame structures, and an- 
holonomically induced torsion even their curvature of the canonical distin- 
guished connection is zero. 

The utihty of the method of anholonomic frames with associated nonhn- 
ear connection structures is that we can work with respect to such frames 
and correspondingly canonical metrics as defined on horizontal/ vertical pairs 
of symmetric Riemannian spaces whose metric and linear connection struc- 
ture is soldered from an underlying Klein geometry. In this regard, the 
bi-Hamiltonian solitonic hierarchies are generated as nonholonomic distri- 
butions of horizontal and vertical moving equations and conservation laws, 
containing all information for a regular Lagrangian and corresponding Euler- 
Lagrange equations. 

We note that curve-flow solitonic hierarchies can be constructed in a sim- 
ilar manner, for instance, for Einstein-Yang-Mills-Dirac equations, derived 
following the anholonomic frame method, in noncommutative generalizations 
of gravity and geometry and possible quantum models based on nonholo- 
nomic Lagrange-Fedosov manifolds. During a long term review of this pa- 
per, a series of new important curve-flow solitonic results were obtained. 
For instance, the results of [22] were generalized in [12] for arbitrary Einstein 
and (pseudo) Riemannian metrics for which alternative d-connections (to 
the Levi Civita one) were defined. It was proven that such metric compat- 
ible d-connections with constant curvature and Ricci d-tensor coefficients, 
with respect to certain N-adapted basis, can be constructed. All construc- 
tions can be redefined equivalently for the Levi Civita connection (because 
all connections under consideration are completely defined by a metric struc- 
ture). 

Finally, we conclude that solitonic hierarchies and the bi-Hamilon and N- 
connection formalisms happen to very efficient in the theory of nonholonomic 
Ricci flows and evolution of physically valuable nonlinear gravitational wave 
solutions [^ mi Us] and Fedosov quantization of Einstein and Lagrange- 
Finsler spaces [16], HT]. So, the results of this paper can be naturally general- 
ized for arbitrary classical and quantum (commutative and noncommutative 
of type [23]) nonholonomic gravitational interactions, generalized Lagrange- 
Finsler systems and their Ricci flow evolutions (see also review PH]). In 
general, we can encode the information on such field/mechanical/evolution 
models into corresponding solitonic hierarchies and, inversely, to extract cer- 
tain nonlinear interaction and/or evolution models from systems of solitonic 
equations. 
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A Some Local Formulas 

This Appendix outlines some local results from geometry of nonlinear con- 
nections (see Refs. |2l E], [251 [26] for proofs and details). There are two types 
of preferred linear connections uniquely determined by a generic off-diagonal 
metric structure with n + m splitting, see g = g (Bn h iHM : 

1. The Levi Civita connection V = {T/j^} is by definition torsionless, 
T = 0, and satisfies the metric compatibility condition, Vg = 0. 

2. The canonical d-connection F'^^^ = (^L^j^, L'^^, C*^, is also metric 

compatible, i.e. Dg = 0, but the torsion vanishes only on h- and 
v-subspaces, i.e. Tj^^. = and T^"^ = 0, for certain nontrivial values 

of Tj", Tj^j. For simplicity, we omit hats on symbols and write L*^ 

instead of L*^, Tj„ instead of Tj^ and so on, but preserve the general 

symbols D and T'^ap- 

By a straightforward calculus with respect to N-adapted frames and 
([5]), one can verify that the requested properties for D on E are satisfied if 

L]k = ^9''' {Gk9jr + ejQkr - er9jk) , (69) 

LI = ek{N^) + h'''{ekh,-hdcehN^-hdbe,Nt), 

<^jc = ^9'''ec9jk, Cl^ = ^h""^ {eMd + e^Kd - edhc) ■ 

For E = TM, the canonical d-connection D = (/i-D, vD) can be defined in 
a torsionless form0 with the coefficients F^g^ = (^^^j^, Ll^^ 

L'jk = ^9'^{(ik9jh + ejgkh- eh9jk), (70) 
The curvature of a d-connection D, 

7^"^ = DF"^ = dT"^ - r^^ A (7i) 

^^Namely where the d-connection has the same coefhcients as the Levi Civita connection 
with respect to N-elongated bases (jH) and (O 
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pa 

^ bjk 


= ^kL\ 


pi 

jka 


= ^-L^k 


qi 

^ jbc 


= e^C^f, 


qa 
bed 





splits into six types of N-adapted components with respect to (jlj) and 
T?° F?" P* P"^ ^ 

-•^ /3^S y-'^ hjki ^ bjki hjai ^ bjai ^ jbei ^ bdc) i 

' ^jL\k + ^"hj^^mk ~ ^"hk^^mj ~ ^^ha^'^kji {^'^) 
•^J-*^ 6fc "T -'^ bi-*^ cfc ^ bk-^ ej be'''' kji 

CfoO + O j^Ly — O ^-^O 

I /^a /^e /^a 

be ~ ^e^ bd + ^ be^ ed ~ bd^ ec" 

Contracting respectively the components, Rq,/? = R^Q,/3r' o'^^ computes 
the h- v-components of the Ricci d-tensor (there are four N-adapted com- 
ponents) 

Rij = R^ijki Ria ^ ~P\kai Rai ^ P^aib^ ^ab ^ "S"^^^^. (73) 

The scalar curvature is defined by contracting the Ricci d-tensor with the 
inverse metric g"^, 

R = g°^R,^ = g'm,, + h-'Sab =R + 'S. (74) 
If E =TM, there are only three classes of d-curvatures, 

^j^\k + ^"hj^\k ~ L^hk^^'mj ~ ^\a^"'kj^ {^^) 

DfeC^ja + ^''jb'P\a^ 
/^a /^a I /^e /^a /^e /^a 

be ^e^ 6d "T 6c'-' ed fed'-' eel 

where all indices a,b, . . . ,i, j, . . . run over the same values and, for instance, 
—>■ C^jk, etc. 
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